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In this article, we study the mass spectrum of the scalar hidden charm and 
Q_i! bottom tetraquark states with the QCD sum rules. The numerical results 

' are compared with the corresponding ones from a relativistic quark model 

■^■j- ■ based on a quasipotential approach in QCD. The relevant values from the 

CN . constituent diquark model based on the constituent diquark masses and the 

' spin-spin interactions are also discussed. 

PACS number: 12.39.Mk, 12.38.Lg 
Q-i| Key words: Tetraquark state, QCD sum rules 

^ : 1 Introduction 

> : 

The Babar, Belle, CLEO, DO, CDF and FOCUS collaborations have discovered (or 
confirmed) a large number of charmonium-like states X(3940), X(3872), F(4260), 
F(4008), F(3940), F(4325), F(4360), ^(4660), etc, and revitalized the interest 
in the spectroscopy of the charmonium states [U El El S] • F° r a concise review 
of the experimental situation of the new charmonium-like states, one can consult 
Ref. [5]. Many possible assignments for those states have been suggested, such as 
J> \ multiquark states (irrespective of the molecule type and the diquark-antidiquark 

type), hybrid states, charmonium states modified by nearby thresholds, threshold 
cusps, etc [TJ [2j [3j 0]. The observed decay channels are J/ipn + 7r~ or tp'n + n~, an 
essential ingredient for understanding the structures of those mesons is whether or 
not the 7T7T comes from a resonance state. 

The Z + (4430) observed in the decay mode ifj'7t + by the Belle collaboration is 
the most interesting subject [6]. We can distinguish the multiquark states from the 
hybrids or charmonia with the criterion of non-zero charge. The Z + (4430) can't 
be a pure cc state due to the positive charge, and may be a ccud tetraquark state. 
The Babar collaboration did not confirm this resonance [7J, i.e. they observed no 
significant evidence for a Z(4430) signal for any of the processes investigated, neither 
in the total J / ipir or ip'n mass distribution nor in the corresponding distributions for 
the regions of Kir mass for which observation of the Z(4430) signal was reported. 

In 2008, the Belle collaboration reported the first observation of two resonance- 
like structures (thereafter we will denote them as Z(4050) and Z(4250) respectively) 
in the tt + Xci invariant mass distribution near 4.1 GeV in the exclusive decays B° — > 
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K tt + Xci [§]• Their quark contents must be some special combinations of the ccud, 
just like the Z + (4430), they can't be the conventional mesons. They may be the 
tetraquark states P [TO] or the molecular states dH E21 E31 EU ■ The Z(4050) and 
Z(4250) lie about (0.5 — 0.6) GeV above the n + Xci threshold, the decay Z — > n + Xd 
can take place with the "fall-apart" mechanism and it is OZI super-allowed, which 
can take into account the large total width naturally. 

The spins of the Z(4050) and Z(4250) are not determined yet, they can be 
scalar or vector mesons. In Refs.JHl [TO] , we assume that the hidden charm mesons 
Z(4050) and Z(4250) are vector (and scalar) tetraquark states, and study their 
masses with the QCD sum rules. The numerical results indicate that the mass 
of the vector hidden charm tetraquark state is about Mz = (5.12 ± 0.15) GeV or 
Mz = (5. 16 ±0.16) GeV, while the mass of the scalar hidden charm tetraquark state 
is about Mz = (4.36 ± 0.18) GeV. The scalar hidden charm tetraquark states may 
have smaller masses than the corresponding vector states. 

The mass is a fundamental parameter in describing a hadron, whether or not 
there exist those hidden charm tetraquark configurations is of great importance 
itself, because it provides a new opportunity for a deeper understanding of the low 
energy QCD. 

In this article, we study the mass spectrum of the scalar hidden charm and 
bottom tetraquark states using the QCD sum rules [TjjJ [16]. In the QCD sum 
rules, the operator product expansion is used to expand the time-ordered currents 
into a series of quark and gluon condensates which parameterize the long distance 
properties of the QCD vacuum. Based on the quark-hadron duality, we can obtain 
copious information about the hadronic parameters at the phenomenological side 

nans]. 

The Z(4050) and Z(4250) can be tentatively identified as the scalar hidden charm 
(cc) tetraquark states, while the scalar hidden bottom (bb) tetraquark states may be 
observed at the LHCb, where the bb pairs will be copiously produced. The hidden 
charm and bottom tetraquark states (Z) have the symbolic quark structures: 

Z+ = QQud; Z° = -^-=QQ(uu - dd); Z~ = QQdu; 

v2 

%t = QQ u s'i Z~ = QQsu; Z® = QQds; ~Z? S = QQsd; 

Zp = —=QQ{uu + dd); Z^ = QQss , (1) 
v2 

where the Q denotes the heavy quarks c and b. 

The colored objects (diquarks) in a confining potential can result in a copious 
spectrum, there maybe exist a series of orbital angular momentum excitations. In 
the heavy quark limit, the c (and b) quark can be taken as a static well potential, 
which binds the light quark q to form a diquark in the color antitriplet channel. 
We take the diquarks as the basic constituents following Jaffe and Wilczek [TTf [TBI. 
The heavy tetraquark system could be described by a double-well potential with 
two light quarks q'q lying in the two wells respectively. 
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The diquarks have five Dirac tensor structures, scalar C75, pseudoscalar C, vec- 
tor 6*7^75, axial vector C*7 M and tensor Ca^ u . The structures Oy^ and Ca^ v are 
symmetric, the structures C75, C and 6*7^75 are antisymmetric. The attractive in- 
teractions of one-gluon exchange favor formation of the diquarks in color antitriplet 
3 C , flavor antitriplet 3/ and spin singlet l s [HI [20]. In this article, we assume the 
scalar hidden charm and bottom mesons Z consist of the C75 — C75 type diquark 
structures rather than the C — C type diquark structures, and construct the inter- 
polating currents: 
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(2) 



where the i, j, k, ■ ■ ■ are color indexes. In the isospin limit, the interpolating currents 
result in three distinct expressions for the correlation functions H(p) , which are 
characterized by the number of the s quark they contain. 

The article is arranged as follows: we derive the QCD sum rules for the scalar 
hidden charm and bottom tetraquark states Z in section 2; in section 3, numerical 
results and discussions; section 4 is reserved for conclusion. 



2 QCD sum rules for the scalar tetraquark states 

Z 

In the following, we write down the two-point correlation functions H(p) in the QCD 
sum rules, 

n(p) = ijd 4 xe ip - x (0\T{j(x)f(0)}\0), (3) 

where the J(x) denotes the interpolating currents Jz+(x), Jz°(x), J z +(x), etc. 

We can insert a complete set of intermediate hadronic states with the same 
quantum numbers as the current operator J(x) into the correlation functions H(p) 
to obtain the hadronic representation [T5l [16]. After isolating the ground state 
contribution from the pole term of the Z, we get the following result, 
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where the pole residue (or coupling) \ z is defined by 

\ z = (0\J(0)\Z(p)) . (5) 

After performing the standard procedure of the QCD sum rules, we obtain the 
following six sum rules: 

\e m 1 = / dspi(s)e , (6) 
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where the i denote the ccqq, ccqs, cess, bbqq, bbqs and bbss channels, respectively; the 
s° are the corresponding continuum threshold parameters and the M 2 is the Borel 



parameter; a max 
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rhq = -^jZa)- The thresholds Aj can be sorted into three sets, we introduce the qq, 
qs and ss to denote the light quark constituents in the scalar tetraquark states to 
simplify the notations, A qq - = ArriQ, A qS = (2rriQ + m s ) 2 , A sS = 4(m<2 + m s ) 2 . 

We carry out the operator product expansion to the vacuum condensates adding 
up to dimension-10. In calculation, we take assumption of vacuum saturation for 
high dimension vacuum condensates, they are always factorized to lower condensates 
with vacuum saturation in the QCD sum rules, factorization works well in large N c 
limit. In this article, we take into account the contributions from the quark conden- 
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sates, mixed condensates, and neglect the contributions from the gluon condensate. 
The contributions from the gluon condensates are suppressed by large denominators 
and would not play any significant roles for the light tetraquark states [211 E2] , the 
heavy tetraquark state [10] and the heavy molecular state [23]. There are many 
terms involving the gluon condensate for the heavy tetraquark states and heavy 
molecular states in the operator product expansion (one can consult Refs. [lOl [23] 
for example), we neglect the gluon condensates for simplicity. Furthermore, we ne- 
glect the terms proportional to the m u and m^, their contributions are of minor 
importance. 

Differentiate the Eq.(6) with respect to jp, then eliminate the pole residues Aj, 
we can obtain the six sum rules for the masses of the Z, 

J A \ dspi{s)e 



3 Numerical results and discussions 

The input parameters are taken to be the standard values (qq) = — (0.24±0.01 GeV) 3 , 
(ss) = (0.8 ± 0.2) (qq), (qg s aGq) = m 2 (qq), (sg s aGs) = m 2 (ss), m 2 = (0.8 ± 
0.2) GeV 2 , m s = (0.14±0.01) GeV, m c = (1.35±0.10) GeV and m b = (4.8±0.1) GeV 
at the energy scale about \i = 1 GeV [T5l [161 121] ■ 

In the conventional QCD sum rules [151 US]; there are two criteria (pole domi- 
nance and convergence of the operator product expansion) for choosing the Borel pa- 
rameter M 2 and threshold parameter sq. The light tetraquark states can not satisfy 
the two criteria, although it is not an indication non-existence of the light tetraquark 
states (For detailed discussions about this subject, one can consult Refs.[I~0~l 125]). 
We impose the two criteria on the heavy tetraquark states to choose the Borel 
parameter M 2 and threshold parameter s - 

If the resonance-like structures Z(4050) and Z(A250) observed by the Belle col- 
laboration in the n + Xci invariant mass distribution near 4.1 GeV in the exclusive 
decays B° — ► K~n + Xci are scalar tetraquark states [8] , the threshold parameter can 
be tentatively taken as s° 5 = (4.248 + 0.5) 2 GeV 2 w 23 GeV 2 to take into account 
all possible contributions from the ground states, where we choose the energy gap 
between the ground states and the first radial excited states to be 0.5 GeV. Taking 
into account the SU (3) symmetry of the light flavor quarks, we expect the threshold 
parameters and s% are slightly larger than the s^. Furthermore, we take into 
account the mass difference between the c and b quarks, the threshold parameters in 
the hidden bottom channels are tentatively taken as = 138 GeV 2 , s° qS = 140 GeV 2 
and s° sS = 142 GeV 2 . 

Here we take it for granted that the energy gap between the ground states and 
the first radial excited states is about 0.5 GeV, and use those values as a guide to 
determine the threshold parameters so with the QCD sum rules. 
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The contributions from the high dimension vacuum condensates in the operator 
product expansion are shown in Figs. 1-2, where (and thereafter) we use the (qq) 
to denote the quark condensates (qq), (ss) and the (qg s aGq) to denote the mixed 
condensates (qg s aGq), (sg s aGs). From the figures, we can see that the contributions 
from the high dimension condensates change quickly with variation of the Borel 
parameter at the values M 2 < 2.6 GeV 2 and M 2 < 7.0 GeV 2 for the cc channels and 
bb channels respectively, such an unstable behavior can not lead to stable sum rules, 
our numerical results confirm this conjecture. At the values M 2 > 2.6 GeV 2 and 
So > 23 GeV 2 , the contributions from the (qq) 2 + (qq) (qg s aGq) term are less than (or 
equal) 10% for the ccqq channel, the corresponding contributions are smaller for the 
ccqs and cess channels; the contributions from the vacuum condensate of the highest 
dimension (qg s ^Gq) 2 are less than (or equal) 2% for all the cc channels, we expect 
the operator product expansion is convergent in the cc channels. At the values 
M 2 > 7.0 GeV 2 and sq > 136 GeV 2 , the contributions from the (qq) 2 + (qq)(qg s (?Gq) 
term are less than 10% for the bbqq channel, the corresponding contributions are 
smaller for the bbqs and bbss channels; the contributions from the vacuum condensate 
of the highest dimension (qg s &Gq) 2 are less than (or equal) 6% for all the bb channels, 
we expect the operator product expansion is convergent in the bb channels. 

In Fig. 3, we plot the contributions from different terms in the operator product 
expansion. From the figures, we can see that the main contributions come from the 
perturbative term and the (qq) + (qg s aGq) term, the operator product expansion 
is convergent; and the interpolating currents contain more s quarks have better 
convergent behavior. 

In this article, we take the uniform Borel parameter M^ in , i.e. M^ in > 2.6 GeV 2 
and M^ in > 7.0 GeV 2 for the cc channels and bb channels, respectively. 

In Fig.4, we show the contributions from the pole terms with variation of the 
Borel parameter and the threshold parameter. The pole contributions are larger than 
(or equal) 50% at the value M 2 < 3.2 GeV 2 and s > 23 GeV 2 , 23 GeV 2 , 24 GeV 2 
for the ccqq, ccqs, cess channels respectively, and larger than (or equal) 50% at the 
value M 2 < 8.0 GeV 2 and s > 136 GeV 2 , 138 GeV 2 , 138 GeV 2 for the bbqq, bbqs 
and bbss channels respectively. Again we take the uniform Borel parameter M^ ax , 
i.e. M^ ax < 3.2 GeV 2 and M^ ax < 8.0 GeV 2 for the cc channels and bb channels, 
respectively. 

In this article, the threshold parameters are taken as Sq = (24 ± 1) GeV 2 , (24 ± 
1) GeV 2 , (25 ± 1) GeV 2 , (138 ± 2) GeV 2 , (140 ± 2) GeV 2 and (140 ± 2) GeV 2 for the 
ccqq, ccqs, cess, bbqq, bbqs and bbss channels, respectively; the Borel parameters are 
taken as M 2 = (2.6 - 3.2) GeV 2 and (7.0 - 8.0) GeV 2 for the cc channels and bb 
channels, respectively. In those regions, the two criteria of the QCD sum rules are 
full satisfied PHES]. 

Taking into account all uncertainties of the input parameters, finally we obtain 
the values of the masses and pole resides of the Z, which are shown in Figs. 5-6 and 
Tables 1-2. 

From Table 1, we can see that the SU(3) breaking effects for the masses of 
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Figure 1: The contributions from different terms with variation of the Borel param- 
eter M 2 in the operator product expansion. The A and B denote the contributions 
from the (qq) 2 + {qq){qg s &Gq) term and the (qg s aGq) 2 term, respectively. The (I), 
(II) and (III) denote the ccqq, acqs and cess channels, respectively. The notations 
a, P, 7, A, p and r correspond to the threshold parameters so = 21 GeV 2 , 22 GeV 2 , 
23 GeV 2 , 24 GeV 2 , 25 GeV 2 and 26 GeV 2 , respectively. 
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Figure 2: The contributions from different terms with variation of the Borel param- 
eter M 2 in the operator product expansion. The A and B denote the contributions 
from the {qq) 2 + (qq) (qg s aGq) term and the (qg s aGq) 2 term, respectively. The (I), 
(II) and (III) denote the bbqq, bbqs and bbss channels, respectively. The notations a, 
(3, 7, A, p and r correspond to the threshold parameters sq = 132 GeV 2 , 134 GeV 2 , 
136 GeV 2 , 138 GeV 2 , 140 GeV 2 and 142 GeV 2 , respectively. 
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Figure 3: The contributions from different terms with variation of the Borel pa- 
rameter M 2 in the operator product expansion. The A, B, C, D, E and F de- 
note the ccqq, ccqs, cess, bbqq, bbqs and bbss channels, respectively. The a, (3 
and 7 correspond to the perturbative term, the (qq) + (qg s crGq) term and the 
(qq) 2 + (qq) (qgscGq) + ('iQsO'Gq) 2 term, respectively. The threshold parameters 
are Sq = 24 GeV 2 and 138 GeV 2 for the cc channels and bb channels, respectively. 
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Figure 4: The contributions from the pole terms with variation of the Borel pa- 
rameter M 2 . The A, B, C, D, E and F denote the ccqq, ccqs, cess, bbqq, bbqs and 
bbss channels, respectively. In the cc channels, the notations a, /3, 7, A, p and r 
correspond to the threshold parameters s = 21 GeV 2 , 22 GeV 2 , 23 GeV 2 , 24 GeV 2 , 
25 GeV 2 and 26 GeV 2 respectively ; while in the bb channels they correspond to the 
threshold parameters s = 132 GeV 2 , 134 GeV 2 , 136 GeV 2 , 138 GeV 2 , 140 GeV 2 and 
142 GeV 2 respectively. 
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Figure 5: The masses of the scalar tetraquark states with variation of the Borel 
parameter M 2 . The A, B, C, D, E and F denote the ccqq, ccqs, cess, bbqq, bbqs 
and bbss channels, respectively. 
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tetraquark states 


This work 


Refs.[26, 27J 


cess 


4.44 ±0.16 


4.051 


acqs 


4.39 ±0.16 


3.922 


ccqq 


4.37 ±0.18 


3.812 


bbss 


11.31 ±0.16 


10.662 


bbqs 


11.33 ±0.16 


10.572 


bbqq 


11.27 ±0.20 


10.471 



Table 1: The masses (in unit of GeV) for the scalar tetraquark states. 



tetraquark states 


pole residues (lO^GeV 5 ) 


cess 


4.00 ±0.80 


ccqs 


3.56 ±0.70 


ccqq 


3.41 ±0.70 


bbss 


17.3 ±4.0 


bbqs 


17.4 ±4.0 


bbqq 


15.2 ±3.8 



Table 2: The pole residues for the scalar tetraquark states. 

the hidden charm and bottom tetraquark states are buried in the uncertainties. 
The central value of the scalar tetraquark state ccqq is slightly larger than the one 
Mz = (4.36±0.18) GeV obtained in Ref.[10j. where the contributions from the terms 
involving the gluon condensate are taken into account. We can draw the conclusion 
that the gluon condensate plays a tiny important role and can be safely neglected. 

The meson Z(4250) may be a scalar tetraquark state (ccud), the decay Z(4250) — > 
7r+ Xci can take place with the OZI super-allowed "fall-apart" mechanism, which 
can take into account the large total width naturally. Other possibilities, such as a 
hadro-charmonium resonance and a D^D° + D + D\ molecular state are not excluded; 
more experimental data are still needed to identify it. It is difficult to identify the 
Z(4050) as the scalar tetraquark state {ccud) considering its small mass. There still 
lack experiential candidates to identify the scalar tetraquark states ccqs, cess, bbqq, 
bbqs and bbss. 

In Table 1, we also present the results from a relativistic quark model based 
on a quasipotential approach in QCD [26], [27], the central values of our predic- 
tions are larger than the corresponding ones from the quasipotential model about 
(0.4 — 0.7) GeV. In Refs. [26| [27] . Ebert et al take the diquarks as bound states 
of the light and heavy quarks in the color antitriplet channel, and calculate their 
mass spectrum using a Schrodinger type equation, then take the masses of the di- 
quarks as the basic input parameters, and study the mass spectrum of the heavy 
tetraquark states as bound states of the diquark-antidiquark system. In the con- 
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ventional quark models, the constituent quark masses are taken as the basic input 
parameters, and fitted to reproduce the mass spectra of the well known mesons and 
baryons. However, the present experimental knowledge about the phenomenological 
hadronic spectral densities of the tetraquark states is rather vague, whether or not 
there exist tetraquark states is not confirmed with confidence, and no knowledge 
about the high resonances. The predicted constituent diquark masses can not be 
confronted with the experimental data. 

In Refs. [251 [251 [3D], Maiani et al take the diquarks as the basic constituents, 
examine the rich spectrum of the diquark-antidiquark states with the constituent 
diquark masses and the spin-spin interactions, and try to accommodate some of 
the newly observed charmonium-like resonances not fitting a pure cc assignment. 
The predictions depend heavily on the assumption that the light scalar mesons 
a (980) and / (980) are tetraquark states, the basic parameters (constituent diquark 
masses) are estimated thereafter. The predications M c5q q = 3723 MeV [28] and 
M C css = 3834 MeV [30] (for the tetraquark states ccqq and cess respectively) are 
about 0.6 GeV smaller than the corresponding ones in the present work. 

In Ref . [3Tj . Zouzou et al solve the four-body (QQqq) problem by three different 
variational methods with a non-relativistic potential considering explicitly virtual 
meson-meson components in the wave-functions, search for possible bound states 
below the threshold for the spontaneous dissociation into two mesons, and observe 
that the exotic bound states QQqq maybe exist for unequal quark masses (the ratio 
mQ/m q is large enough). The studies using a potential derived from the MIT bag 
model in the Born-Oppenheimer approximation support this observation [32, 33J. 
In Ref. [33], Manohar and Wise study systems of two heavy-light mesons interact- 
ing through an one-pion exchange potential determined by the heavy meson chiral 
perturbation theory and observe the long range potential maybe sufficiently attrac- 
tive to produce a weakly bound two-meson state in the case Q = b. In Ref.[35j, 
the L = tetraquark states QQQQ (Q denotes both Q and q) are analyzed in a 
chromo-magnetic model where only a constant hyperfine potential is retained. 

If there exist scalar tetraquark states QQqq, we can construct the C-y^ — Cj 1 * 
type interpolating currents to study them with the QCD sum rules, as the QQ and 
qq favor forming diquarks in the symmetric sextet 6/ with the spin-parity J p = 1 + 
due to Fermi statistics. The attractive interactions of one-gluon exchange favor 
formation of the diquarks in color antitriplet 3 C , flavor antitriplet 3/ and spin singlet 
l s [19l[20]. We expect the scalar Cj^ — Cj 11 type tetraquark states QQqq are heavier 
than the corresponding C75 — C75 type tetraquark states QQqq, our numerical 
results support this conjecture; our works on the C*7 M — Cj^ type tetraquark states 
QQqq will be presented elsewhere. 

The nonet scalar mesons below 1 GeV (the / (980) and a (980) especially) are 
good candidates for the tetraquark states. However, they can't satisfy the two 
criteria of the QCD sum rules, and result in a reasonable Borel window. If the 
perturbative terms have the main contribution (in the conventional QCD sum rules, 
the perturbative terms always have the main contribution), we can approximate 
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the spectral density with the perturbative term [25], then take the pole dominance 
condition, and obtain the approximate relation, 

£>4.7. (11) 

If we take the Borel parameter has the typical value M 2 = lGeV 2 , then so > 
4.7GeV 2 , the threshold parameter is too large for the light tetraquark state candi- 
dates / (980), a (980), etc. 

On the other hand, the numerous candidates with the same quantum numbers 
J PC = ++ below 2 GeV can't be accommodated in one qq nonet, some are supposed 
to be glueballs, molecules and multiquark states [HJ EE1 [37] . Once the main Fock 
sates of the nonet scalar mesons below 1 GeV 2 are proved to be tetraquark states, 
we can draw the conclusion that the QCD sum rules are not applicable for the light 
tetraquark states. 

In this article, we calculate the mass spectrum of the scalar hidden charm and 
bottom tetraquark states by imposing the two criteria of the QCD sum rules. In fact, 
we usually consult the experimental data in choosing the Borel parameter M 2 and 
the threshold parameter sq. There lack experimental data for the phenomenological 
hadronic spectral densities of the tetraquark states, the present predictions can't be 
confronted with the experimental data. 

The LHCb is a dedicated b and c-physics precision experiment at the LHC (large 
hadron collider). The LHC will be the world's most copious source of the b hadrons, 
and a complete spectrum of the b hadrons will be available through gluon fusion. 
In proton-proton collisions at i/i = 14 TeV, the bb cross section is expected to be ~ 
producing 10 12 bb pairs in a standard year of running at the LHCb operational 
luminosity of 2 x 10 32 cm _2 sec _1 [3H]- The scalar tetraquark states predicted in the 
present work may be observed at the LHCb, if they exist indeed. We can search for 
the scalar hidden charm tetraquark states in the DD, D*D*, D S D S , D*D*, J/ipp, 
J/ip4>, J/i/juj, i] c ir, rj c 7], ■ ■ ■ invariant mass distributions and search for the scalar 
hidden bottom tetraquark states in the BB, B*B*, B S B S , B*B*, Tp, T0, Tu, rj^Tr, 
?7{,77, • • • invariant mass distributions. 

Furthermore, the non-leptonic B decays through b — *> ccs provide another favor- 
able environment for the production of the scalar hidden charm tetraquark states 
[39], we can search for them at the KEK-B or the Fermi-lab Tevatron. 



4 Conclusion 

In this article, we study the mass spectrum of the scalar hidden charm and bottom 
tetraquark states with the QCD sum rules. The numerical results are compared with 
the corresponding ones from a relativistic quark model based on a quasipotential 
approach in QCD. The relevant values from the constituent diquark model based on 
the constituent diquark masses and the spin-spin interactions are also discussed. We 
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can search for the scalar hidden charm and bottom tetraquark states at the LHCb, 
the KEK-B or the Fermi-lab Tevatron. 
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